Semi-Inclusive A and Kg Production in p-Au Collisions at 17.5 GeV/c 
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The first detailed measurements of the centrality dependence of strangeness production in p-A 
collisions are presented. A and Kg dn/dy distributions from 17.5 GeV/c p-Au collisions are shown as 
a function of "grey" track multiplicity and the estimated number of collisions, v, made by the proton. 
The v dependence of the A yield deviates from a scaling of p-p data by the number of participants, 
increasing faster than this scaling for v < 5 and saturating for larger v. A slower growth in Ks 
multiplicity with v is observed, consistent with a weaker v dependence of KK production than YK 
production. 
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Significant effort has been devoted in the last decade 
to measuring strange particle production in nucleus- 
nucleus (A- A) collisions, motivated by predictions that 
quark-gluon plasma (QGP) formation could enhance 
strangeness Experiments at both the Brookhaven 

National Laboratory (BNL) AGS and the CERN SPS 
accelerators have reported large increases in relative 
strange particle yields in central light (Si, S) f^-|6| and 
heavy ion (Au, Pb) J7|— pTT]| induced collisions compared 
to p-p collisions. However, we still cannot claim a true 
understanding of the physics of strangeness enhancement 
due to the complexity of the hadronic interactions un- 
derlying A-A collisions and the competing mechanisms 
proposed and/or used in models Jl2|-p^| to explain the 
data. 

The difficulty in interpreting the A-A data and the 
observation that the enhancement is already present in 
light-ion collisions suggest the use of p-A collisions to 
study this problem further. The simpler final state of 
p-A collisions may allow the production rate for strange 
particles to be directly connected with the scattering dy- 
namics of the incoming proton. Previously published p-A 
data Jl5|,|l6| demonstrated an increase in the inclusive 
K/ir ratio with increasing A at the AGS, suggesting that 
the strangeness enhancement mechanism is already at 
work in p-A collisions. A more thorough analysis JLj} of 
inclusive data at higher energies |l^] suggested no overall 
strangeness enhancement but a possible modest enhance- 
ment in A production offset by a decrease in K produc- 
tion. This result has led to claims that the observed 
strangeness enhancement in A-A collisions at SPS ener- 
gies may result from QGP formation |19| . However, ex- 
trapolations to central A-A data from inclusive p-A data 
are intrinsically flawed: Centrally selected A-A events 
necessarily involve more scatterings of the participant nu- 
cleons, and the dynamics of strangeness production may 
be quite sensitive to these additional scatters. A resolu- 
tion of this problem requires a detailed measurement of 
the centrality dependence of strangeness production in 
p-A collisions. 

In this paper we present the first such measurement, 
made by BNL experiment 910 at the AGS accelerator. 
The data consist of A and Kg rapidity spectra and in- 
tegrated yields obtained as a function of "grey" track 
multiplicity from p-Au collisions at a beam momentum 
of 17.5 GeV/c. E910's nearly complete rapidity coverage 
allows us to accurately estimate the total A and Kg mul- 
tiplicities and study the variation of absolute yields with 
v(N grey ), the estimated number of collisions of the pro- 
ton in the Au nucleus. A common benchmark fl ^flll 
for evaluating strangeness enhancement in A-A collisions 
is the scaling of p-p data by the number of participants, 
./Vpart, in the collision. For p-A collisions, an -ZVp art scaling 
of p-p data would yield [£lj 

N P ro d = \N^ od (l + u), (1) 

since there are two participants in p-p collisions. We com- 



pare our multiplicities to the yields expected from Eq. 1 
to evaluate whether we see enhanced strange particle pro- 
duction using the same benchmark as in A-A collisions. 
We observe that for A production Eq. 1 is physically sen- 
sible since baryon number is conserved and BB processes 
are negligible at our energies. Deviations from Eq.l 
would imply contributions from target nucleons not di- 
rectly struck by the projectile or changes in the probabil- 
ity of participants to fragment into A's. E910 was staged 
in the Multi-Particle Spectrometer (MPS) facility at the 
AGS. For this data, the Cherenkov tagged secondary pro- 
ton beam had a mean momentum of 17.5 ± 0.1 GeV/c 
and a 1.5% momentum spread. The E910 spectrometer 
was previously described in ^,^3| ; the results presented 
here rely on the EOS time projection chamber and the 
beamline/trigger counters. A 3.9 g/cm 2 Au target was lo- 
cated 20.5 cm upstream of the active area of the TPC and 
was immediately followed by a 10 cm x 10 cm two-layer 
scintillating-fiber hodoscope. This multiplicity detector 
provided two triggers, a minimum-bias (MB) trigger that 
required 2 hits on each layer and a "central" trigger that 
required a total of 20 hits in the two layers and selected 
approximately the 25% most central events. Due to the 
low light yield of the fibers the MB trigger suffered sig- 
nificant efficiency loss for low multiplicity events with no 
highly ionizing tracks. 

After finding and fitting the recorded pulses in the 
TPC we obtained typical resolutions of 0.7 mm (verti- 
cal) and 0.5 mm (horizontal) for position measurements 
in each sample. The momentum resolution for particles 
bending in the 0.5 T magnetic field varied from 1.2% (p < 
2 GeV/c) to 5.5% (p ~ 17 GeV/c). The dE/dx measure- 
ment, obtained from a truncated-mean of the TPC sam- 
ples on each track, provided a resolution of a /(dE/dx) 
= 6% for typical track lengths. A's and Kg's were mea- 
sured and identified through a combination of topolog- 
ical reconstruction and dE/dx identification of the de- 
cay daughters. We paired and removed conversion elec- 
trons and positrons, the dominant source of background, 
with an efficiency of « 50%. We further reduced back- 
ground from conversions and from false vertices by apply- 
ing tighter geometric cuts to small opening-angle pairs. 
We attempted topological fits on +/— track pairs satis- 
fying applied geometric and dE/dx cuts and accepted as 
Vo candidates those passing applied y 2 cuts with origin 

> 3.5 cm from the target. We used a combined like- 
lihood from the single track dE/dx's and hypothetical 
p?r~ and 7r + 7r~ invariant masses (Mj nv ) of the pairs to 
identify the decaying particle, obtaining the Mj nv distri- 
butions shown in Fig. [y. We obtained mass resolutions of 
FWHM=4 MeV, 13 MeV for A's and Kg's, respectively 
and corresponding S/B ratios of 35:1 and 30:1. The ac- 
ceptances were calculated via GEANT simulations of the 
detector response to 20M pure A and 10M pure Kg de- 
cays [p4| . We show in Fig. [I] the y — p± regions with 

> 10%acceptance. For the centrality measurement, we 
identified as "grey" tracks protons and dcuterons in the 
momentum ranges [0.25,1.2] GeV/c and [0.5,2.4] GeV/c, 



respectively. We obtain the multiplicity of grey tracks, 
iVgrgy, within our geometric acceptance event-by-event 
and from this quantity estimate the mean number of col- 
lisions suffered by the beam proton, v(N gTey ), using an 
established technique p3| . 

The data presented in this paper resulted from a com- 
bined 4.65M triggers. We required a valid event to have 
at least one secondary charged particle in the final state 
and a T,p± > 85 MeV/c. In addition, we vetoed one 
and two-track events containing a high-momentum pos- 
itive track consistent with a quasi-elastically scattered 
proton. After applying quality and the above interac- 
tion cuts we obtained 2.97M events, 1.88M minimum- 
bias and 2.07M central. From these, we reconstructed a 
total of 156.8k A's and 76.8k Ks's. Using beam-triggered 
events we determined trigger efficiency corrections as a 
two-dimensional function of charged-particle multiplicity 
and Agroy . The correction for multiplicity 1,2 events with 
Agrcy = is large (5.1) while the average correction for 
all interactions is 1.1. 

We calculated A and K s yields AN(m±, y, M inv ) per 
event as a function of A grey after subtracting the A and 
Ks background in each bin. We corrected these for ac- 
ceptance (A), trigger efficiency (e), and branching ratio 
(BR) to obtain an invariant differential yield, 

d 2 n A/Ks _ 1 1 AN A/Ks (m ± ,y) 
m±dm±dy m±AeBR, N ovt Am±Ay 

For Ag rcy < 4, we used only MB triggers, but for larger 
A gr cy we combined both triggers, weighting by the num- 
ber of events in each sample. The m± spectra are uni- 
formly well-described by exponential distributions except 
at the highest m± values (toj_ — m > 0.6 at mid-rapidity) 
where we do not have sufficient statistics to accurately 
determine the background subtraction. We fit the m± 
spectra excluding these points to the form, 

1 d2jl = 1 dn c -(m ± -m )/T^ (3) 

2nm± dm±dy 27r(mo + T)T dy 

where dn/dy is a direct parameter of the fit represent- 
ing the integral of Eq. ^| over m± . The fits give inverse 
slopes, T, that vary from 0.05 GeV/c at low and high y to 
0.14 GeV/c at mid-rapidity, consistent with proton spec- 
tra obtained from p-A collisions at a similar energy [p.6| . 
For A grcy bins where we do not have enough statistics 
to perform the fits we directly sum over mi to obtain 
dn/dy. Fig. || shows the resulting dn/dy distributions 
for a sub-set of the available A groy bins. We estimate 
90% CL point-to-point systematic errors in the A and 
Kg dn/dy measurements, including contributions from 
the fitting, to be < 5% except for the lowest rapidity bin 
(20%), and estimate the normalization systematic error 
to be ±10%. 

We observe that with increasing A groy , the A and Kg 
yields decrease at high rapidity and increase at low ra- 
pidity. For A's, the decrease in yield at large rapidity 



is a direct consequence of the increased "stopping" of 
the projectile baryon resulting from the multiple interac- 
tions in the target nucleus ^2j. The strong increase in 
production of strange particles at low rapidity is qualita- 
tively consistent with previously observed trends in sec- 
ondary particle production in p-A collisions fl6| , ^5| , p6f . 
We show in Fig. || the integrated yields that we have ob- 
tained by summing our measured dn/dy values as a func- 
tion of v. The error on the in-acceptance yield is domi- 
nated by and taken to be the same as our uncertainty in 
the absolute normalization ±10% at 90%CL. We show 
(solid line) in Fig. || the expectations from Eq. [j] with 
N£ p = 0.054±0.002(syst) and = 0.035±0.002(sys£) 
obtained by parameterizing the ^/s dependence of A and 
Kg multiplicities |27| and interpolating to our energy. 
The A yields initially increase faster with v than expected 
from the A part scaling of p-p yields and then saturate and 
start to decrease. The Kg yields behave similarly with 
a slower initial increase. The apparent decrease of the 
yields may result from the fact that we miss a larger frac- 
tion of the total A and Kg yield with increasing A grcy or 
v due to our low-rapidity cut-off. We have estimated the 
missing yield by fitting the dn/dy distributions to gamma 
distributions as shown in Fig. |^ and extrapolating these 
into the unmeasured region to produce the estimated to- 
tal yields shown in Fig. |^. The uncertainty in the total 
yield is largest for the larger A grcy or v bins where the 
dn/dy distribution peaks near the edge of our acceptance. 
We show in Fig. || 90%CL systematic errors on the yields 
with larger errors on the high side to account for the 
possibility that an unknown mechanism may produce a 
larger A yield below y = than we estimate. The result- 
ing total A yields shown in Fig. || saturate at large v and 
remain flat. We have fit the extrapolated A yields to an 
empirical function, 

A A =A p A p (l-e— Q )/(l- e ~ K ), (4) 

where N pp is the A multiplicity in p-p collisions. The 
obtained function with n = 0.299 ±0.008 and a = 1.29 ± 
0.03 describes both the initial rapid rise of the yield and 
the saturation at large v. To evaluate the significance of 
this fast initial increase in the A yield we plot in Fig. || the 
yield that would result from a "binary-collision" scaling 
of p-p data, 

N BC {v) = N pp v, (5) 

which we view as the fastest plausible increase that could 
be expected from the multiple scattering of the incoming 
proton. The A yields are consistent with this "upper 
limit" for v < 3 indicating a rate of increase in A yield 
with v in this region that is approximately twice that 
given by Eq. 1 . We note that the systematic error on the 
v scale of ±15% [^3| is small compared to this difference 
in slope. 

We observe that the Kg yield increases more slowly 
with v than the A yield and appears to decrease slightly 
at large v even after we have accounted for the missing 



yield. The difference in behavior between the A and K5 
yields may result from the mixture of K° and K° in the 
Ks and the fact that these are produced through differ- 
ent processes. In p-p collisions at comparable energies, 
A of the K5 are produced as K°'s in association with 
hyperons with the other half produced as K°'s associ- 
ated with kaons If we assume that this proportion 
is not modified in p-A collisions and that the total hy- 
peron yields increase in proportion to the A yield, we can 
estimate the K° component of the measured K5 yields 
shown in Fig. |^. The K° component appears to increase 
by a factor of two at v = 3, roughly consistent with the 
-Wpart scaling of p-p data shown in Fig. ||, before starting 
to decrease slowly with v. While our data suggests an 
increase in KK production with v, because of uncertain- 
ties in the above assumptions we cannot make a stronger 
statement. A forthcoming analysis of K~ production will 
provide clearer insight into this problem. 

In conclusion, we have reported on the first detailed 
investigation of the centrality dependence of strange par- 
ticle production in p-A collisions using grey track mul- 
tiplicity and the estimated number of collisions of the 
projectile nucleon to characterize centrality. We have 
measured dn/dy distributions for A and K5 that both 
show a strong backward shift with increasing N grcy and 
v. The estimated total A yields increase with v at a rate 
approximately twice that expected from the Af part scaling 
of p-p data for v < 3 and saturate for v > 5. As noted 
above, this violation of Af part scaling implies either that 
additional nucleons not directly struck by the projectile 
contribute to A production or that the probability for one 
or more of the participants to fragment into a A increases 
with v. We observe a slower but significant increase in Kg 
multiplicity with v that apparently results from different 
behavior of the K° and K° components of the K5. The 
observed increase in Kg yield is large enough to allow for 
a statistically significant increase in KK production with 
v for v < 3 using a reasonable extrapolation of p-p data. 

We conclude that at AGS energies, p-A data show a 
clear violation of a simple Y part scaling of p-p data. This 
result has clear qualitative implications for use of such 
scaling for interpreting strangeness yields in A-A colli- 
sions. To quantitatively evaluate the potential implica- 
tions of our results, we assume that the target contri- 
bution to the p-Au A yield grows as vN pp /2 and at- 
tribute the remainder to the fragmentation of the pro- 
jectile and/or energy deposition of the projectile in the 
nucleus. Then, our data show that the "projectile" con- 
tribution increases proportional to v for v < 3 with a 
slope that is the same as for the target nucleons. In A-A 
collisions where both the projectile and target nucleons 
multiply scatter, this picture implies that the hyperon 
and associated kaon yields per participant would increase 
rapidly with the average number of scatters of the par- 
ticipants, (v), for (y) < 3 giving a maximum possible 
increase in yield per participant of a factor ~ 3. This is 
precisely the behavior seen in the K + production in Si-Au 



and Au-Au collisions at the AGSJpj and A production 
in Pb-Pb collisions at the SPS JTlJT In p-A collisions, 
the enhancement is more modest - a 50% increase over 
A^part scaling at v = 3 - simply because the target nucle- 
ons scatter only once. As we have shown, however, this 
modest enhancement may have profound consequences 
for interpretation of the strangeness enhancement in nu- 
clear collisions. We note that the above picture picture 
is consistent with the additive quark model [ p9[ used by 
Kadija et al. fl30] | to explain the observed strangeness en- 
hancement in light-ion collisions at the SPS. In partic- 
ular, the increase in the projectile-like component up to 
v = 3 is exactly what is expected from the additive quark 
model. Since the saturation of the A yield for v > 5 is 
very likely due to the stopping of the incident baryon we 
predict that at higher energies the A yield will continue 
to increase for v > 5. 
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FIG. 3. A (left) and K s (right) multiplicities vs u(N srey ). 
T - in E910 acceptance, A - total yield (syst. errors from 
extrapolation included), o - est. K° component of Ks. Lines: 
Solid - iVp a rt scaling, dot-dashed - Binary collision, dashed - 
empirical fit to A (left panel, see text for details), est. -/V part 
scaling for K° (right panel). 
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FIG. 1. A and Ks acceptance boundary (left) and recon- 
structed invariant mass spectra (right). 
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FIG. 2. A (left) and Ks (right) rapidity density distribu- 
tions for different Agrcy plotted in both log (top) and lin- 
ear(bottom) scale. A grcy = • - 0, o - 1, T-2 and V - 4. 



